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ABSTRACT 

A projected mass estimator and a maximum likelihood estimator are developed to 
best determine the mass of stellar clusters from proper motion data. Both methods 
effectively account and correct for errors in the velocities, and provide an unbiased 
and robust estimate of the cluster mass. Using an extensive proper motion study of uj 
Centauri by van Leeuwen et al. (2000), we estimate the mass of u Centauri using these 
two different methods to be 2.8 x 10 ®Mq [d/b.\kpc\ . Within this modeling context, 
the statistical error is 3% and the systematic error is 4%. 
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1 INTRODUCTION 

Estimating accurate masses of self-gravitating systems, such 
as star clusters, is of central importance in Astronomy. Most 
estimates of such masses have been based on radial velocities 
from a few tracers of the system. However, the advent of the 
next generation telescopes, proper motion data will become 
more accurate and reliable within an observation time of a 
few years. Proper motions at the Milky Way’ s center has 
already proven a boon le.g. ISchodel et ^1200,'^) . Combined 
with radial velocities, proper motions can provide us with 
5 components of the 6-dimensional phase space of a large 
number of tracers. Such a wealth of information will help to 
set accurate limits to masses of objects in the local group. 
Apart from the masses, this new kinematic data can uncover 
more subtle dynamic processes in star clusters, helping us 
to verify some of our present understanding and models of 
star clusters. 

The idea of using proper motion data for cluster dynam¬ 
ics is not new. Kudworth was the first to implement this idea 
in the 70’s, but he was limited by the data and telescopes 
of the time. iLeonard fc MerrittI dlQSflll showed that proper 
motion data was sufficient to predict the dynamic mass of 
the cluster. However, there was no complete proper motion 
data sample of a star cluster at the time. 

At the present, data from HST have reached sufficient 
quality to study internal cluster k i nemat ics astrometrically 
[McNamara. Harrison fc Anderson] (l20n.'il . Also, diffraction 
limited imaging at the VLT, Keck and other 8 m class tele- 
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scopes is now a reality, and with the next generation of tele¬ 
scopes c oming online (e.g. the LET with 23 m resolution 
imaging, iHerbst et a O|200i), proper motions are becoming 
a broadly competitive dynamical tool, and are not restricted 
to a reas with long time-base line data, such as lo Centauri 
(e.g. Ivan Leeuwen et alJ^OOflll . Hence, it is crucial to under¬ 
stand how best to use proper motion data for setting mass 
limits and to learn to deal with the systematic and random 
errors inherent to this technique. 

In Section we undertake a review of present mass 
estimators and develop tools geared towards proper motion 
data, in the hope of extracting maximum information about 
the same. We then apply these tools to a proper motion data 
set of u! Centauri (Figure^ in Sectionj^to predict the mass 
of the cluster. 


2 MASS ESTIMATORS 

The virial theorem has been traditionally and widely used 
to estimate the masses of self-gravitational systems from the 
motions of individual tracers or constituents. However, a 
major drawback of estimators that are based on the virial 
theorem is that t hey tend to be biased and fo rmally have an 
infinite variance. iBahcall fc Tremainel (Il98lll compared the 
virial mass estimator to an alternative method called the 
projected mass e stimator first proposed bylPagj tl95^ for 
binary galaxies. iHeisler. Tremaine fc Bahcffi J |I198 !tI 1 later 
expanded their work to develop the projected mass estima¬ 
tor for self-gravitating systems. 

While for many self-gravitating systems (e.g. extra- 
galactic ones) radial velocities are the only measurable quan- 
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Figure 1. Observational c onstraints for proper motion based mass estimates of lo Centauri : (a) Logarithm of the observed surface 
brightness profile taken from iMevlanI il987t) . (b) Radial (solid circles) a nd transverse (open triangles) proper motion dispersion profiles as 
a function of projected radius, as derived bv ivan Leeuwen et alJ j20QOli from the complete proper motion data set. The central dispersion 
profile is uncertain due to crowding effects at the centre of the cluster (indicated by the shaded region), (c) The errors in the radial 
(solid circles) and transv erse (open triangles) proper motions of the chosen restricted sample of class 0 and class 1 stars (4828 stars. See 
Ivan Leeuwen et al ]|20 o 3) used to determine the mass of the cluster. 


titles, it has become increasingly clear in the last decade that 
proper motions stu dies could be used t o estim ate the masses 
of nearby clusters. lT,eonard fc MerritH ill QSijl extended this 
approach to include proper motions in the mass estimate of 
the open cluster M35. 

2.1 Projected Mass Estimator (PME) 

Projected mass estimators try to exploit the directly ob¬ 
servable, projected properties of a self-gravitating system. 
In addition, they are efficient, unbiased and have a hnite 
variance. The estimators are derived by taking moments of 
the Jeans equation for a spherical system. 

The projected mass estimator for proper moti on veloc¬ 
ities can be rewritten as (iLeonard fc MerrittlllQSiil Eg. 19). 

M = — <{2vl+vl)R> (1) 

where vr and vt are the radial and transverse proper motion 
velocities, R is the projected distance and M is the total 
mass of the cluster. 

Since the projected mass estimator involves a line-of- 
sight projection, care must be taken to have a proper radial 
sampling of the tracers of the system. Note that in the pro¬ 
jected mass estimator, particles at larger projected distances 
are given more weight than particles at smaller projected 
distances, because of the inability to reconstruct the true 
radial distance for particles at small projected distance. 

The projected mass estimator becomes biased when 
the velocity measurements have signihcant error bars. The 
estimated mass bias is a quadratic effect w.r.t the error 
in velocity. However, we can modify the projected mass 
estimator to take into account the velocity errors. Since 
< ^obs >=< Vtrue > + < Rtrror >: ^ reasonable unbiased 
estimator is 

Op 

^ ~ 3tvG ^ +Vt — Uer-ror)^ >i (2) 

where v'^rror = ^Svji + SVr. 


The projected mass estimator does not provide a “for¬ 
mula” for mass error estimates. However bootstrapping 
methods can be used to set limits to the mass. 

Correcting for the errors in the velocities as in the above 
estimator is a statistical process, and may fail with a lim¬ 
ited number of particles. Overestimation or underestimation 
of the errors in the proper motion velocities may bias the 
estimator. 


2.2 Maximum Likelihood Estimator (MLE) 

Here we develop a maximum likelihood method mass estima¬ 
tor that is geared towards proper motion measurements and 
accounts for signihcant, and individually varying, error bars. 
It uses high-precision measurements effectively and extracts 
maximum information from low-precision velocity measure¬ 
ments. The key to this method is a suitable model to pre¬ 
dict the proper motion velocity dispersions (both radial and 
transverse) as a function of the projected radius. For exam¬ 
ple, the velocity dispersion can be predicted by the family of 
King models or by Jeans modelling of the cluster. The data 
then simply set the mass normalization of the model. 

Let us consider that we can model the proper motion 
velocity distribution locally as a Maxwellian distribution. 
Thus the probability of a star in the cluster to have a certain 
proper motion velocity v given that the total mass of the 
cluster is M is: 

Pmodei{v\M) = — ^ expi-— ^ 1 (3) 

y/2-K aRM,R) \2aRM,R)j 

where the velocity dispersions vary with radius and scale 
with mass as a{M,R) = '/M<j{R) and R is the projected 
radial distance. 

On the other hand, every observation of proper motion 
velocity Vi has a certain error 5vi associated with it. The 
probability of actually measuring the velocity Vi at a pro¬ 
jected distance Ri is a Gaussian distribution given by : 
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Thus the probability that the cluster has a mass M 
given an observation of proper motion velocity Ui, at a pro¬ 
jected distance Ri with an observational error Svi is, 

Pi{vi\M) = J Pob3{vi\v,Svi)pmodei{v\M)dv. (5) 

The log-likelihood of the cluster to have a mass M given 
N observations is: 

N 

C{yvi\ M) = ^ lnpi{vi\ M), (6) 

i=l 

where one sums over both radial and transverse proper mo¬ 
tion velocities. 

The most probable mass of the cluster corresponds to 
the parameter M of the maximum of the likelihood func¬ 
tion. Limits of uncertainty of Icr can also be set to the mass 
of the cluster by searching for the parameter M such that 
SC = Cmax — T = 0.5 (in the astrophysical context see 


Lampton. Margon & Bowver 

IT97«. whose statistic S = 2C\ 

Efstathiou. Ellis & PetersonI 

1988|1. 


3 APPLICATION TO uj CENTAURI 

The Southern globular cluster u) Centauri (N GC 5139) 

a proximately at a distance of 5.1 kpc dLongmorel 
. Its large mass and size sets it apart from the 
bulk of globular clusters in the Galaxy. In recent years, 
it has evoked interest because of the recent discovery of 
its multiple populations which implies a complex forma- 
tion history I^Uke^^^^cMi^lioOTj^nghe^^^^Wte^ei^ 
|2000j}/rheduster is also isotropicT ilKhi^^^^ideraoDl^nn^L 
iMemt^^^J (|W9^ reported a definite rotation in the clus¬ 
ter. T^ideraM ll2002ll reports that there is very little mass 
segregation in the cluster, which means that the cluster is 
not in equipartition. 

There has also been a considerable debate on the mass 
of UJ Gentauri. Previous estimate s of the mass of uj Genta uri 
have ranged fr om 2.4 x 10^Mp llMandushev et alJll99lll to 
7.13 X 10® Mq llRlcher et ahlioolll . Setting accurate limits to 
the mass of uj Centauri is essential to differentiate between 
the various theories of i t s form ation process. 

Ivan Leeuwen et (l2nnff> have completed a project 
providing the largest available set of proper motion data 
for a star cluster. This large data set on uj Centauri made 
available offers a fore-taste of the proper motion data which 
may be soon made available for more clusters. 

3.1 Data set available 

For our mas s estimates, we use the pr oper motion studies of 
UJ Centauri llvan Leeuwen et^ ]|200^, which are based on 
100 plates obtained with the Yale-Columbia 66 cm refractor, 
and ranging in epoch from 1931 till 1935 and 1978 till 1983, 
differential proper motions are obtained for 9847 stars. The 
limiting magnitude (photographic) is 16.0 for the centre of 
the cluster, 16.5 for the outer part. The precisions of the 
proper motions range from an average of 0.1 mas/yr for the 
brightest to an average of 0.65 mas/yr for the faintest stars. 


Ivan Leeuwen et alJ i200fll uses the proper motions to study 
the internal kinematics of the cluster. 

The centre of the cluster was taken to be q = 
201°.69065, 5 = —47°.47855. We assume that the distance of 
the cluster is 5.1 kpc away iLongmorell992^ . We project the 
absolute celestial coordinates a and S into Cartesian coor¬ 
dinates (relative to the cluster centre) onto the plane of the 
sky along the line-of-sight vector through the cluster centre: 

X = —/o cos (5 sin Aa 

y = /o (sin (5 cos 5o — cos (5 sin (5o cos Aa) (7) 

where /o is a scaling factor (/o = 1 to have x and y in radians 
and /o = 180/7r to have x and y in units of degrees), ao and 
So are the cluster centre coordinates, and Aa = a — ao and 
/\S = S-So. _ 

Ivan Leeuwen et al.l (l200(t> divide their proper motion 
observations into 4 classes based on the errors in the proper 
motion determination which depend upon crowding and 
plate quality. They also calculate the membership proba¬ 
bility of each star to belong to the cluster. 

We only use stars taken from high-quality classes 0 and 
1, since t he errors appear to be rela tively well behaved (See 
Fig. 11 of ivan Leeuwen ~ 'all l2nnfii . 

The mass estimators are very sensitive to outliers and 
field stars in the cluster. Thus the mass estimates show a 
direct correlation to our membership probability cut-off cri¬ 
terion. We use a conservative cut-off criterion: we use stars 
having a membership probability higher than or equal to 
99.73% {3a) to belong to the cluster. This gives us a lim¬ 
ited sample of 4056 stars, which have a 100% probability to 
belong to the cluster, due to the integral nature of the mem - 
bership probability assigned bv Ivan Leeuwen et alJ i200Cll . 
These we term as class A stars. 

There are 772 stars having a membership probability 
between 98.5% {2.5a) and 99.73% to belong to the cluster. 
These are class B stars. To diminish the effect of any non¬ 
members on the mass estimators, we construct a Monte- 
Carlo test by using any subset of class B stars along with 
class A stars. In choosing our subset, we ensure that the 
number of class B stars does not exceed 5% of the class A 
stars. 

As shown in Figure Q the errors in the proper motion 
data becomes comparable to the velocity dispersion profile 
at a magnitude of 16, especially for stars in the outer parts 
of the cluster. Thus reasonable and robust estimates of the 
mass can be obtained at magnitude of 15. 

Recent proper motion s observations of uj Centauri 
by iKing fc AndersonI ^2^^2^ with the WFPC2 aboard the 
HST have shown that the data from Ivan Leeuwen et alJ 
i2QQ0l) may be p r oblem atic at the centre of the cluster. 
iKing fc AndersoM i2002ll obtain proper motions dispersions 
much lower than reporte d in the proper motion study of 
Ivan Leeuwen et alJ 11200011 : however, the two groups agree 
with each other farther out. This will not affect our mass 
estimates from the PME, as it gives a lower precedence to 
data at the centre of the cluster. For the MLE, we neglect 
stars from the inner region of 1/ this is indicated by the 
shaded region of Figure Q 

The radial stellar profile was taken from iMevlanI il987ll . 
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Figure 2. Results for the Projected Mass Estimator (Eq. O as a 
function of the class 0 and 1 stars from ivan Leeuwen et 
The 1 — a error bars are calculated using bootstrapping. The 
closed symbols use class A stars only. The open symbols use both 
class A and class B stars in a Monte-Carlo method. 


3.2 Projected Mass Estimator 

In this section, we attempt to estimate the mass of uj Cen- 
tauri f rom the proper motion data set of ivan Leeuwen et alJ 
i2000ll using the tools we have developed in the previous 
section. First, we apply the PME (Eq. to the data set; 
then we use the MLE from Section f2.21 with various models 
of the projected velocity dispersion profile. 

One must be wary of selection effects in applying the 
PME. Crowding effects and lack of colour information at the 
centre of the cluster leads to a bias in the radial sampling of 
stars with proper motions. A simple check can be made by 
calculating the average projected radius as a function of the 
limiting magnitude for our restricted sample; we find that 
fainter sub-samples lie at larger radii. Thus we have a larger 
number of faint stars away from the centre of the cluster in 
the sample. 

We can correct for this effect by using a Monte-Carlo 
algorithm to choose our tracers for the PME such that we 
always maintain the surface density profile of the cluster. 
For this purpose, we fit a Kings densit y profile to sur face 
density data (see Figure 0 published bv iMevlanI il98i1l . 

The 1-cr error bars (see Figure 121 are calculated us¬ 
ing bootstrapping. At very bright magnitudes, we have high 
precision measurements, but only few stars. The large error 
bars demonstrate the lack of information due to the limited 
amount of tracers of the system. 

In trying to extract the maximum information in an 
economic way with the PME, one has to trade off between 
the above two effects. At magnitude of 15 using only class A 
stars, the PME estimates the mass of the cluster as 2.72 x 
10® Mq, while the combination of class A as well as class B 
stars in a Monte-Carlo method predicts 2.78 x 10 ®Mq. 


3.3 Maximum Likelihood Estimator 

To apply the MLE to oj Centauri, we require a suitable 
model for the v elocity disp ersion. We start with the clas¬ 
sic King model (lKinelll96(ill framework for this task. Given 
a suitable dynamical model of the cluster, the tracers of 
the system effectively scale the projected velocity dispersion 
curve with the help of a maximum likelihood method to es¬ 
timate the independent parameter: the mass of the cluster. 

As with the PME, we also correct for “proper radial” 
sampling. We also neglect all stars in the central region of 
1' due to the uncertainty in the central velocity dispersion 
from crowding effects. At most 1.3% of all the tracers of our 
limited sample are present in this central region. 


3.3.1 King Models 

In Figure 121 a single mass King Model is fit to the surface 
density profile of lo Centauri. The central scaled potential 
used is VFo = 6, to match the surface radial profile llMevlanI 
1198711 . Since the search fo r anisotropy in u) Centa uri has 
yielded very limited results llKine fc Andersorl2002^ . we are 
justified in using one integral models. Even though the sur¬ 
face density profile provides a relatively good fit, we are 
unable to reproduce the projected velocity dispersion pro¬ 
file in the inner parts of the cluster (even well outside !')• 
The resulting mass estimate from the MLE as a function of 
limiting magnitude is also shown. 

At magnitude of 15 using only class A stars, the MLE 
with the King model estimates the mass of the cluster 
as 2.83 X 10® Mq, while the combination of class A as 
well as class B stars in a Monte-Carlo method predicts 
2.96 X 10® Mq. 


3.3.2 Jeans modelling 


Without assuming any restrictive form of a distribution 
function, we can still reconstruct the velocity dispersion us¬ 
ing Jeans equation. Under the hypothesis of isotropy and 
neglecting the rotation in the cluster, the isotropic Jeans 
equation for a spherical system is: 




( 8 ) 


where a is the radial velocity dispersion, v is the observable 
number density distribution, and is the potential gen¬ 
erated by the matter distribution at the considered point. 
Note that uj Centauri is one of the few systems known to be 
isotropic from direct observations. 

W e use the three dimen sional number density distribu¬ 
tion bv iMerritt et al.l (ll997^ using non-parametric analysis. 
Assuming that mass does not follow light, we try to disen¬ 
tangle the mass profile from the luminosity p rofile, by giving 
it an independent King profile ||King|[l2S^- Note that the 
King mass profile only defines the potential well of the clus¬ 
ter, and does not assume the particular form of the Kings 
distribution function. The mass profile was adjusted until 
the form of the velocity dispersion curve predicted by Jeans 
equatio n matched the kinematic d ata of u) Centauri as de¬ 
rived bv lva.n Leeuwen et alJ feOPfil . This is shown in Figure 
0] The velocity dispersion curve is then used to predict the 
mass of the cluster using the MLE. 
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Figu re 3. (a) The Surface brightness profile of the single mass King model (Wo = 6) compared with the sur face density data tM exlaiJ 
Il987^ (b) Corresponding projected velocity dispersion profile of the King model as compared with the data ivan Leeuwen er’aHI^OOQlb 
The dashed vertical line indicates the half-mass radius of this model, (c) The resulting mass estimate of the MLE as a function of limiting 
magnitude. The closed symbols use class A stars only. The open symbols use both class A and class B stars in a Monte-Carlo method. 


The shape of the mass profile is consistent with the fact 
that uj Centauri exhibits very litt le eqnipartition, as shown 
by the results of lAndersonI (|2QQ^, where it was found that 
stars of a considerable range of mass down the main seqnence 
of u! Centauri differed very little in their radial distribution; 
the differences found indicate only a small evolution towards 
eqnipartition. The lack of equipartition is to be expected 
with a cluster of such large mass, and hence a very long 
relaxation time. 

At lower magnitude the error-bars are large due to lack 
of information in the limited number of tracers. At magni¬ 
tude of 15 using only class A stars, the MLE with the Jeans 
model estimates the mass of the cluster as 2.72 x 10®M©, 
while the combination of class A as well as class B stars in 
a Monte-Carlo method predicts 2.85 x 10 ®Mq. 

The resulting MLE mass estimates for the cluster from 
the two models differ by 4% because of the difference in 
the velocity dispersion profiles at the central region of the 
cluster (< 5'). On the other hand, the results of the MLE 
with Jeans modelling are consistent with the results of the 
PME. 


4 CONCLUSIONS 

In this paper, we implement two independent methods to de¬ 
termine the mass of a cluster from the proper motion data. 
When using the PME, one must be careful to correct for the 
errors in the velocities (Eq. and to use a “proper radial 
sampling” to correct for the incompleteness at the centre of 
the cluster. The PME provides us an unbiased and robust 
estimate of the mass of the cluster, so long as the errors 
in the velocities are as large as the velocity dispersion it¬ 
self. Using stars above this limiting magnitude leads us to 
a systematic error. We independently estimated the mass 
of the cluster using the MLE in conjunction with two dy¬ 
namical models for the velocity dispersion profile. In this 
paper, we use the King model and Jeans modelling. Dis¬ 
entangling the mass profile from the density profile in the 
Jeans modelling, gives us a better prediction of the shape 
of the velocity dispersion profile. Both models give us very 


similar mass estimates of the cluster within a 4% systematic 
error. The mass estimators are susceptible to a systematic 
bias due to contamination by outliers and field stars in the 
cluster. 

The largest uncertainty in the absolute mass estimate is 
the distance to u) Centauri which enters as the third power. 
A com parison of the radial velocity dispersion llMavor et alJ 
Il997^ with the proper motion velocity dis persions indicates a 
dista nce of the cluster of 4.74T0. 16 kpc ifvan Lee uwen et ^ 
l200(fl . However we follow Iron Leeuwen et in interpreting 
these results with caution, and assume that the cluster is at 
a distance of 5.1 kpc away. 

Thus, when implemented correctly, the two indepen¬ 
dent methods p redict the mass of the clus ter within 6% on 
the basis of the Ivan Leeuwen et al.l i200flll data. The mass 
of u] Centauri is 2.8 x 10 ®Mq [d/b.lkpd^ with a 3% statis¬ 
tical error and a 4% systematic error which arises from the 
various models used. The choice of the limiting magnitude 
introduces a further systematic error of 4%. 
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solid line. The dashed line is the assumed King mass profile (not normalisedy_£b]_The_corres£ondingvelocity dispersion curve predicted 
by Jeans equation is plotted along with the kinematic data as derived by |va£_Jjeeuwenmt_^^ J200C|). The vertical dotted line on the 
velocity dispersion plot indicates the half mass radius, (c) The resulting mass estimate of the MLE as a function of limiting magnitude. 
The closed symbols use class A stars only. The open symbols use both class A and class B stars in a Monte-Carlo method. 
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